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Ceria (Ce0 2 ) is a promising catalyst for the reduction of carbon dioxide (C0 2 ) to 
liquid fuels and commodity chemicals, in part because of its high oxygen storage 
capacity, yet the fundamentals of C0 2 adsorption and initial activation on Ce0 2 sur- 
faces remain largely unknown. We use density functional theory, corrected for onsite 
Coulombic interactions (DFT+U), to explore various adsorption sites and configura- 
tions for C0 2 on stoichiometric and reduced Ce0 2 (110). Our model of reduced Ce0 2 
(110) contains oxygen vacancies at the topmost atomic layer and undergoes surface 
reconstruction upon introduction of these vacancies. We find that C0 2 adsorption on 
reduced Ce0 2 (110) is thermodynamically favored over the corresponding adsorption 
on stoichiometric Ce0 2 (110). The most stable adsorption configuration consists of 
C0 2 adsorbed parallel to the reduced Ce0 2 (110) surface, with the molecule situated 
near the site of the oxygen vacancy. Structural changes in the C0 2 molecule are 
also observed upon adsorption, so that the resulting O-C-0 angle is 136.9° and the 
C-0 bonds are 1.198 A and 1.311 A in length, respectively. The molecule bends 
out of plane to form a unidentate carbonate, as opposed to the bidentate carbonate 
found by other researchers for CO adsorption to stoichiometric Ce0 2 . We deduce 
that charge transfer from reduced surface Ce 3+ ions to the adsorbate to form the 
carbonate anion is the first step in the activation and reduction process and cleavage 
of the C-0 bond. 
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I. INTRODUCTION 



The synthesis of liquid fuels and commodity chemicals by carbon dioxide (CO2) reforming 
is a promising approach for clean energy production. For this reason, much academic and 
industrial effort has been devoted to exploring efficient means of reducing COj^. It is widely 
believed that the first step in CO2 reduction is the activation of the C=0 bond and eventual 
formation of the CO2 radicaP. Since CO2 is thermodynamically stable and the formation 
of the radical is difficult in the gas phase^ (i.e., -1.90 V vs. NHEP), we seek appropriate 
catalysts for lowering the activation barrier to (photo)electrochemical CO2 reduction. 

There have been several metal oxides proposed for CO2 activation, such as ZrOj^, TiOj^, 
MgQi33 anc i Q a oDJ] Among these various catalysts for CO2 activation, ceria (Ce02) is 
attracting increased interest due to its high oxygen storage capacities. Ce0 2 is commonly 
used for oxidation and reduction reactions, since Ce ably and reversibly converts between 
Ce 4+ and Ce 3+ upon release and storage of oxygerP. 

While industrial methods to efficiently utilize carbon monoxide (CO), which has similar 
chemistry to CO2, have been developed, CO2 has been comparatively underutilized as a 
reactant. We already know that the interaction between ceria and CO is vitally important 
in applications ranging from three-way automotive catalysts (TWC™, the water-gas shift 
(WGS) reaction™ 1 , and syngas production! 16 * 17 *. However, all of the research to date has yet 
to produce a comprehensive model of how ceria and CO2 interact at a molecular leve l 18 * 19 *, 
although there have been several experimental studies with support for metal 

nanoparticles that have hinted at important factors; researchers have noted the importance 
of the interaction between the metal nanoparticle and the oxide support for CO2 activation^ 
and the possibility of oxygen cycling between the support and the adsorbed CO2 species^^ll. 
In fact, CO is formed by the decomposition of CO2 on partially reduced Pt/ceria catalysts, 
whereas on fully pre-oxidized (or stoichiometric) ceria, little to no reactivity is observed^. 
However, these experimental studies did not clarify the mechanism of CO2 activation on 
ceria nor the structure of the adsorbed species on the surface. Thus, it is important to 
investigate CO2 activation on ceria at a molecular level. 

A number of first-principles studies on ceria have explored the stoichiometric bulk and as- 
sociated surfaces. Reduced ceria, with oxygen vacancies at the surface, has also been studied 
extensively^^. It has been reported that the ceria (110) surface is more catalytically active 
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than either the (111) or the (100) surfaces, due to the metastable state of (110) surface, so 
that the ordering of surface stability is (111) I (HO) l (lOO^EHES^ However, the energy of 
formation for a single surface oxygen vacancy (with release of 5O2) on stoichiometric ceria 
(111) is +2.65 eV^, while it is +2.11 eV for the ideal ceria (110) surface^. Although the 
ceria (110) surface is the second most stable surface among the low index ceria surfaces, it 
is thus chemically more active than ceria (111). The creation of oxygen vacancies on the 
ceria (110) surface requires the least amount of energy, which facilitates the formation of a 
reduced surface^ESl. 

In this paper we present a comprehensive study on the adsorption of CO2 on the stoi- 
chiometric and reduced ceria (110) surfaces. A variety of possible adsorption configurations 
on these two surfaces are presented in terms of structures and energetics, along with charge 
transfer and electronic property analyses. Our results provide theoretical evidence that 
the initial adsorption of ceria surfaces induces CO2 activation to ultimately form the anion 
radical CO^, and the adsorption process is sensitive to surface structure and composition. 

II. COMPUTATIONAL METHODS 

All calculations are performed within the framework of density functional theory (DFT)^^, 
as implemented in the Vienna Ab Initio Simulation Package (VASP 5.2)'^^' and the general- 
ized gradient approximation of Perdew, Burke, and ErnzerhoP^ to represent the exchange- 
correlation energy. The Projector- Augmented Wave (PAW) methocP^, with a 400 eV 
energy cutoff, was used to describe the wavefunctions of the atomic cores. The tetrahedron 
method with Blochl corrections ^ was used to set the partial occupancies for the orbitals. 
While several fc-point mesh sizes (e.g. 4 x 4 x 4 up to 11 x 11 x 11) were considered, we 
ultimately used the 6x6x6 T-centered Monkhorst-Pack fc-point mesh for bulk ceria to give 
results that were sufficiently converged (within 1 x 10~ 5 eV using the conjugate gradient 
method). 

The bulk ceria unit cell contains 8 oxygen atoms and 4 cerium atoms in a cubic unit 
cell with the Fm3m symmetry space group; the initial lattice parameter is |a| = 5.411 A^. 
After geometry optimization, the bulk lattice parameter increased slightly to |a| = 5.413 A; 
this compares favorably to previous studies that also follow the DFT approach®'. 

The electronic band structure of bulk ceria contains two energy gaps, as determined by 
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XPS spectra: 1. A 3 eV band gap between the O 2p valence band and the Ce 4/ conduction 
band, and 2. A 6 eV band gap between the valence band and the Ce 5d band^M3 Since it is 
difficult to accurately represent the 4/ states in ceria using conventional DFT methods 25133331 , 
we have chosen to employ the Hubbard U parameter within the GGA+U approach! 50 * 51 *. The 
Hubbard U term acts as an on-site Coulombic interaction to properly localize the electrons 
in these states, which is especially important for the calculations on reduced ceria, with its 
partially occupied 4/ states. We have chosen for our ceria system U e ff = 5 eV, which is 
consistent with recommended values previously published in the literature^. 

Slab models consisting of 10 atomic layers and 15 A of vacuum between periodic images 
in the z direction were cleaved along the (110) plane from the optimized bulk CeO"2 unit cell; 
following optimization, we observe that atomic layer relaxation occurred primarily within 
the three topmost layers. Thus, we cleaved the slab in half along the xy plane, so that 
a p (2 x 2) slab consisting of 5 atomic layers (overall slab thickness is 7.622 A), with the 
bottom two layers constrained in space, was used for all subsequent calculations. The other 
atoms are free to move in all directions within the fixed lattice. The exposed surface area 
is large enough to accommodate a small gas-phase adsorbate (e.g., CO2) without spurious 
interactions from neighboring images. 



The optimized stoichiometric CeO"2 (HO) surface is shown in Figure la Several distinct 
adsorption sites are labeled, including: 1. Threefold coordinated O atoms in the topmost 
atomic layer (nO), 2. Fourfold coordinated O atoms in the subsurface atomic layer (nnO), 3. 
Sixfold coordinated Ce atoms in the topmost atomic layer (nCe), 4. Eightfold coordinated 
Ce atoms in the subsurface atomic layer (nnCe), and 5. Hollow. This nomenclature will 
be used throughout this discussion to describe changes in the adsorption configuration and 
electronic structure at the surface. 

The optimized reduced Ce02 (HO) surface is formed by removing one oxygen atom from 
the topmost atomic layer of the p (2 x 2) slab, so that a 12.5% rate of oxygen vacancies is 
created; this is shown in Figure [Tb] Several distinct adsorption sites are labeled, including: 
1. Threefold coordinated O atoms in the topmost atomic layer (nO), 2. Fourfold coordinated 
O atoms in the subsurface atomic layer (nnO), 3. Fivefold and sixfold coordinated Ce atoms 
in the topmost atomic layer (nCe), 4. Sevenfold and eightfold coordinated Ce atoms in the 
subsurface atomic layer (nnCe), 5. Hollow, and 6. Vacancy. We observe that the nO nearest 
the vacancy site moved to fill the vacancy so that it is equidistant between two nCe atoms; 




(a) Stoichiometric (b) Reduced 



FIG. 1: Schematic representations of the CeO"2 (HO) surfaces. The top figures view the 
surfaces from the top, and the bottom figures view the surfaces from the side. nO (large 
red sphere) refers to a surface O atom. nCe (large beige sphere) refers to a surface Ce 
atom. nnO (small red sphere) refers to a subsurface O atom. nnCe (small beige sphere) 

refers to a subsurface Ce atom. 

it also moved out of the surface in the perpendicular direction. Furthermore, the other nnO 
and nnCe atoms near the vacancy moved slightly away (0.07-0.09 A) from the vacancy, so 
the vacancy site is quite large in area and volume. As a result, the nCe-0 distances near 
the vacancy site increase to 2.354 A, compared to 2.273 A in stoichiometric CeO"2 (HO). 
This surface reconstruction in reduced Ce02 (HO) is in agreement with results previously 
reported in the literature®^. 

The energy of adsorption for carbon dioxide on ceria is given as: 

E a ds = E C 2 +Ce0 2 (lW) ~ -^CeO 2 (110) — -E'cCb (1) 

where -Ece0 2 (iio) is the total energy of the Ce02 (HO) surface slab, Eco 2 is the total energy 
of the CO2 molecule (optimized in a periodic cubic unit cell whose volume is 8000 A 3 ), 
and -Eco 2 +CcO 2 (ii0) is the total energy of the composite system. Since the calculations are 
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performed at K and fixed cell volume, the differences in Gibbs free energy should equal 
the differences in total energy. By this definition, a negative value of E ads corresponds to 
an exothermic and spontaneous adsorption process. 

Both parallel and perpendicular /vertical orientations of CO2 on the stoichiometric and 
reduced Ce02 (HO) surfaces were considered, and the energy of adsorption calculated for all 
distinct configurations. Excess spin density and Bader charge^ analyses were also performed 
to study the role of charge transfer in the CO2 adsorption and activation processes at the 
Ce0 2 (HO) surface. 

III. RESULTS AND DISCUSSION 

A. CO2 Adsorption on Stoichiometric CeC>2 (HO) 

The optimized adsorption configurations of CO2 on stoichiometric Ce02 (HO) are shown 
in Figure [2j Seven configurations where CO2 is adsorbed parallel to the surface (along either 
the x- or y-axis) and four configurations where CO2 is adsorbed perpendicular/vertical to 
the surface were considered. The energies of adsorption, select geometrical parameters, and 
variation of the total Bader charge of the adsorbed CO2 molecule are shown in Table |T| 

The strongest energy of adsorption (Eads = —0.262 eV) is observed for nnO adsorbing as 
a "bridge" between two nCe atoms (Sp5). The two C-0 bonds in CO2 are slightly elongated 
to 1.196 A and 1.207 A, compared to 1.176 A in the gas phase. The CO2 molecule are nearly 
linear, with a O-C-0 angle of 178.2°. The two oxygen atoms in CO2 reside at the top of 
two sixfold coordinated surface nCe atoms, and the Ce-0 distances are 2.520 A and 2.527 
A, respectively. An additional depiction of this adsorption configuration (Sp5) is shown in 
Figure |3j 

The Sv4 and Svl adsorption configurations are slightly less energetically favorable 
(E ads = —0.253 eV and E ads = —0.201 eV, respectively), and both are in vertical orien- 
tation and tilted slightly relative to the surface normal. For Sv4, the CO2 adsorbs at the 
nnO_nnO hollow site, which is flanked by two subsurface ceria atoms (nnCe) and two subsur- 
face oxygen atoms (nnO). The distance between the ceria surface and the adsorbing oxygen 
atom in CO2 is 2.631 A. For Svl, the CO2 adsorbs at the nO nO hollow site, which is 
flanked by two surface ceria atoms (nCe) and two surface oxygen atoms (nO). The distance 
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Spl. nO atop-x Sp2. nCe atop Sp3. nnO-nnO bridge 




FIG. 2: CO2 adsorption on p (2 x 2) supercell of stoichiometric Ce02 (HO). Sp denotes 
adsorption parallel to the surface, and Sv denotes adsorption perpendicular /vertical to the 
surface. For the CO2 adsorbate, oxygen is represented by an orange sphere and carbon is 

represented by a black sphere. 

between the ceria surface and the adsorbing oxygen atom in CO2 is 2.573 A. The symmetry 
of these two hollow sites reduces the net intermolecular forces acting on the CO2 molecule, 
so these adsorption configurations are relatively stable. However, the lengths of the two C-0 
bonds in CO2 for the Sv4 and Svl configurations are similar to those in gas-phase CO2. 
Additional depictions of these these two adsorption configurations are shown in Figure [3] 

Several other adsorption configurations merit discussion, in terms of alignment of the 
adsorbate and consideration of surface/subsurface effects. For Spl and Sp4, the C atom 
of CO2 resides on top of the surface nO atom, but is aligned along different axes (x and y, 
respectively). Both Spl and Sp4 result in very weak physisorption, with Sp4 being slightly 
favored energetically due to the increased electrostatic repulsion with the neighboring surface 
nO oxygen in the Spl configuration compared to the repulsion with the subsurface nnO 
oxygens in the Sp4 configuration. For Sp2 and Sv2, where CO2 is adsorbed on top of 
the surface nCe atom, the energy of adsorption is stronger than that for Sp6 and Sv3, 
respectively, where C0 2 is adsorbed on top of the subsurface nnCe atom. Any differences 
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TABLE I: Representative structural parameters, energies of adsorption, and variation of 
the total Bader charge for a CO2 molecule adsorbed at various sites (Figure [2]) on the 

stoichiometric Ce0 2 (HO) surface 



rc-o (A) 


^o-c-o (°) 


Eads (eV) 


Ap(C0 2 ) 

M 


Spl 1.177-1.178 


180 


-0.034 


-0.003 


Sp2 1.182-1.185 


178.8 


-0.118 


-0.002 


Sp3 1.183-1.186 


179.2 


-0.183 


-0.005 


Sp4 1.178-1.182 


180 


-0.078 


-0.005 


Sp5 1.196-1.207 


178.2 


-0.262 


-0.008 


Sp6 1.180-1.182 


179.5 


-0.043 


-0.002 


Sp7 1.179-1.183 


178.1 


-0.174 


-0.002 


Svl 1.172-1.180 


178.2 


-0.201 


-0.006 


Sv2 1.173-1.178 


179.1 


-0.102 


-0.003 


Sv3 1.175-1.176 


180 


-0.035 


-0.001 


Sv4 1.172-1.179 


180 


-0.253 


-0.003 



in electrostatic repulsion with subsurface nnCe or nnO atoms appear to be minimal, as 
the energies of adsorption for the Sp3 and Sp7 configurations are nearly identical. Thus, 
we propose that only surface atoms appear to significantly influence the strength of CO2 
adsorption to the stoichiometric CeC>2 (HO) surface. 

In general, all of these adsorption configurations on stoichiometric CeC>2 (HO) would be 
classified as weak physisorption, since the C0 2 molecule remains physically unchanged upon 
adsorption and the distance between C0 2 and the surface is more than 2.52 A, which is 
too far to form a Ce-0 bond. Furthermore, the Bader charge analysis indicates that the 
variation of the total charge of adsorbed CO2 is less than —0.01 |e|, so there is no significant 
change in the electronic structure. We can conclude that stoichiometric CeC>2 (HO) does 
not activate the C0 2 molecule. 
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Sp5 Sv4 Svl 

f= -0.262 eV £= -0.253 eV E= -0.201 eV 

FIG. 3: CO2 adsorption on stoichiometric Ce02 (HO) in selected configurations: Sp5, 
Sv4, and Svl. Both structural parameters and energies of adsorption are depicted. 



B. C0 2 Adsorption on Reduced Ce0 2 (110) 

The optimized adsorption configurations of CO2 on reduced CeO"2 (HO) are shown in 
Figure |ij Seven configurations where CO2 is adsorbed parallel to the surface (along either 
the x- or y-axis) and four configurations where CO2 is adsorbed perpendicular/vertical to 
the surface were considered. The energies of adsorption, select geometrical parameters, and 
variation of the total Bader charge of the adsorbed CO2 molecule are shown in Table [XTJ 

The strongest energy of adsorption (E ads = —1.233 eV) is observed for C0 2 parallel 
adsorption as a "bridge" between two nCe atoms (Rp5) adjacent to the vacancy site. From 
Figure [4j we can see the carbon atom of the adsorbate is situated directly above the nO atom 
that had moved toward the vacancy oxygen atom. Furthermore, the adsorbed CO2 molecule 
exhibits a bent configuration, where the O-C-0 angle is 136.9°. Also, the two C-0 bonds 
in the adsorbed CO2 molecule are elongated to 1.246 A and 1.256 A, compared to 1.176 A 
in the gas phase. Thus, we would classify this process as chemisorption. Although there 
is significant lengthening of the C-0 bonds, complete dissociation to form CO+O is not 
observed. In addition, the C-0 distance between adsorbed CO2 and the surface nO atom 
is 1.498 A, so a carbonate group appears likely but does not actually form upon surface 
adsorption. An additional depiction of this adsorption configuration (Rp5) is shown in 
Figure |5j 

The next strongest energy of adsorption (E a d s = —1.012 eV) is observed for CO2 par- 
allel adsorption on top of the vacancy (Rp4). The two C-0 bonds in the adsorbed C0 2 
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FIG. 4: C0 2 adsorption on p (2 x 2) supercell of reduced Ce0 2 (HO). Rp denotes 
adsorption parallel to the surface, and Rv denotes adsorption perpendicular /vertical to 
the surface. For the CO2 adsorbate, oxygen is represented by an orange sphere and carbon 

is represented by a black sphere. 




Rp5 Rp4 

F= -1.223 eV E= -1.012 eV 

FIG. 5: CO2 adsorption on reduced Ce02 (HO) in selected configurations: Rp5 and Rp4. 
Both structural parameters and energies of adsorption are depicted. 
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TABLE II: Representative structural parameters, energies of adsorption, and variation of 
the total Bader charge for a CO2 molecule adsorbed at various sites (Figure [ij on the 

reduced Ce0 2 (110) surface 



rc-O (A) 


/-o-c-o (°) 


Eads (eV) 


Ap(C0 2 ) 

M 


Rpl 1.176-1.183 


178.3 


-0.401 


-0.273 


Rp2 1.182-1.197 


135.5 


-0.119 


-0.223 


Rp3 1.207-1.287 


177.3 


-0.201 


-0.119 


Rp4 1.198-1.311 


140.4 


-1.012 


-0.936 


Rp5 1.246-1.256 


136.9 


-1.223 


-0.955 


Rp6 1.179-1.183 


180 


-0.055 


-0.016 


Rp7 1.177-1.182 


178.9 


-0.177 


-0.059 


Rvl 1.168-1.182 


177.1 


-0.815 


-0.081 


Rv2 1.201-1.322 


139.9 


-0.265 


-0.921 


Rv3 1.172-1.181 


179.1 


-0.038 


-0.043 


Rv4 1.175-1.182 


178.5 


-0.282 


-0.178 



molecule are 1.198 A and 1.311 A, respectively. The adsorbed CO2 molecule exhibits a 
similar binding geometry and bent configuration with Rp5; however, the O-C-0 angle is 
slightly larger (140.4°) than the corresponding O-C-0 angle in Rp5 (136.9°). Furthermore, 
the C-0 distance between adsorbed CO2 and the vacancy nO atom in Rp4 is also larger 
(1.642 A), than the corresponding C-0 distance in Rp5, which suggests that the adsorbed 
CO2 molecule is not as activated as in the Rp5 configuration. An additional depiction of 
this adsorption configuration (Rp4) is shown in Figure [5j 

Vertical adsorption atop the vacancy site (Rvl) is slightly less energetically favorable 
(E a d s = —0.815 eV) than for Rp5 or Rp4, and we observe no analogous activation of the 
CO2 molecule; the resulting O-C-0 angle is 177.1° and no significant changes are measured 
in the C-0 bond lengths compared to the gas phase values. Though one oxygen atom 
of the CO2 molecule is inserted into the surface vacancy, the reduced surface is sufficiently 
reconstructed so that the remaining surface oxygen has moved to bridge the adjacent surface 
cerium atoms; thus, the CO2 molecule in the Rvl configuration cannot fully "heal" the 
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vacancy to regenerate the stoichiometric surface. 

Vertical adsorption atop the nCe atom nearest the vacancy site (Rv2) is less energetically 
favorable (E a( i s = —0.265 eV) than adsorption atop the vacancy site (Rvl), although the 
adsorbed molecule is bent similarly to Rp5. The other adsorption configurations on reduced 
CeC>2 (HO), where CO2 is adsorbed far from the vacancy site, are less energetically stable, 
with energies of adsorption on the same order of magnitude as those on stoichiometric CeC>2 
(110); thus, we will not explore these configurations further, and instead, primarily focus 
our subsequent discussion and analysis on the most stable adsorption configuration, Rp5. 



C. Electronic analysis of CeC>2 as a catalyst for CO2 activation 

We now comment on the structural changes induced by oxygen vacancy formation at 
reduced Ce02 (HO) surfaces. We see that for the Rp5 configuration, the Ce-0 distances 
between either neighboring nCe atom and the nO atom at the vacancy are significantly 
longer (2.292-2.237 A) than the corresponding Ce-0 distance in the clean reduced surface 
(2.165 A). This suggests that structural changes occur even following CO2 adsorption, and 
may be attributable to charge transfer between the surface and adsorbed CO2. 

The partial electronic density of states (PEDOS) for Ce and O, when we have the clean 



reduced Ce02 (HO) surface, is shown in Figure 6a, and it is consistent with those previously 
published in the literatur e! 32 * 33 * 58 ! We aligned the top of the valence band to eV, and 
observed a filled gap state in the middle of the band gap; this filled gap corresponds to the 
unoccupied Ce 4/ state. The integration of this peak indicates that there are two electrons 
per oxygen vacancy. To locate these two electrons, we plotted the excess spin density, which 
is defined as the difference between the up spin density and the down spin density, for clean 



reduced Ce02 (HO), as shown in Figure 7a We see that the spin density on the clean reduced 
Ce02 (HO) surface is localized on the two nCe atoms on the surface nearest the vacancy. 
The shapes of the two isosurfaces resemble those of / orbitals, which is consistent with the 
4/ electronic configuration in Ce 3+ . Since cerium exists in the Ce 4+ state in stoichiometric 
Ce02, the excess spin density plot suggests that it adopts the Ce 3+ state upon formation 
of the oxygen vacancy. This result is also in agreement with previous observations on other 
reduced ceria surfaces^. 

The PEDOS for Ce and O in the situation where C0 2 is adsorbed to reduced Ce0 2 
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FIG. 6: Partial electronic density of states (PEDOS) for Ce and O in reduced Ce02 (HO). 
The top of the valence band is aligned to eV, and the Fermi level is indicated with a 

vertical dashed line. 



(110) is shown in Figure 6b We observe that the gap between the Fermi level (top of 
the unoccupied Ce 4/ peak) and the bottom of the conduction band widens upon C0 2 
adsorption, and nears the value observed for stoichiometric Ce02 (110)P^. In addition, the 
filled gap state peak narrows upon CO2 adsorption; the integration of this peak indicates 
that there is only one electron per oxygen vacancy, compared to two electrons in the clean 
reduced surface. To locate the remaining electron, we plotted the excess spin density for 
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FIG. 7: Excess spin density in reduced CeO"2 (HO). Upon formation of the oxygen 
vacancy, adjacent cerium atoms are reduced to Ce 3+ . Upon CO2 adsorption, one cerium 
cation is re-oxidized back to Ce 4+ and the excess charge is transferred to the adsorbate. 
The Ce-0 bonds at the surface are also lengthened. 



reduced CeO"2 (HO) with adsorbed CO2 (Rp5), as shown in Figure 7b From Figure 7b 
we observe that the remaining electron is localized on the nCe2 atom; thus one reduced 
Ce 3+ cation (nCe2) remains on the surface, while the other (nCei) has been re-oxidized 
back to Ce 4+ . The asymmetric electric field at the surface vacancy, with one Ce 3+ and one 
Ce 4+ , explains the asymmetric C-0 bond elongation found in the most stable adsorption 
configuration (Rp5). We also observe that the re-oxidized Ce 4+ atom is closer in distance 
to the nO oxygen atom that moved towards the vacancy site upon surface reduction, while 
the Ce 3+ atom is further from this nO oxygen atom and the adsorbed CO2 molecule due to 
electrostatic repulsion. These results indicate that the reduced ceria (110) surface has been 
partially oxidized upon C0 2 adsorption at the vacancy site. 
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Since the CO2 molecule may accept electrons into its lowest unoccupied molecular orbital 
to form the carbonate anion, we performed a Bader charge analysis on the system, focusing 
on the interface between C0 2 and the reduced Ce0 2 (HO) surface, and plotted the electron 
localization function (ELF)PP. From the variation of the Bader charges on the CO2 molecule 
in various adsorption configurations (Table Table [XT]) , we observe definitive charge transfer 
from the ceria surface to the adsorbed CO2 molecule. This is especially prominent in the case 
of the Rp5 configuration shown in Figure [7| where a net charge of —0.955 |e| is transferred 
to form an adsorbed unidentate carbonate species. This is confirmed by the ELF plot shown 
in Figure |8j where it is clear that the adsorbed CO2 is activated with a net negative charge 
localized on the oxygen atoms of the molecule. Thus, C0 2 is activated to form a unidentate 
carbonate anion upon adsorption parallel to vacancy sites on reduced CeO"2 (HO) surfaces. 




FIG. 8: Electron localization function for reduced CeC>2 (HO) with adsorbed CO2 (Rp5). 
A net charge of —0.955 |e| is transferred from the surface to form an adsorbed unidentate 

carbonate species. 

We also performed a Bader charge analysis for the perpendicular/vertical adsorption con- 
figurations of C0 2 to the reduced Ce0 2 (HO) surface. We observed that while no significant 
charge transfer occurs for the Rvl configuration, the variation in the Bader charge for the 
Rv2 configuration is —0.921 |e|, which again suggests that Ce 3+ is oxidized back to Ce 4+ 
when a CO2 molecule available for reduction is placed in close proximity. However, the low 
energy of adsorption for the Rv2 configuration suggests that it is not favored compared to 
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the Rp5 configuration. 

In summary, we conclude that the specific structure of the ceria surface has an critical 
influence on the CO2 activation. CO2 may chemisorb to reduced CeC>2 (HO) with an energy 
of adsorption of -1.223 eV, and that charge transfer from the reduced ceria surface to the 
adsorbate results in the formation of an activated unidentate carbonate species. In the 
CO2 adsorption process, the reduced ceria (110) surface is only partially re-oxidized, which 
suggests that there is additional capacity for further reaction if the CO2 concentration at 
the surface is increased. Our present findings are in agreement with recent experimental 
work, in which the researchers showed that re-oxidation of ceria surfaces using CO2 would 
only be energetically feasible at the adsorption sites that have been sufficiently reduced^. 
Furthermore, partial re-oxidation of reduced ceria occurs with high reaction probability even 
at room temperature^, so we propose that reduced ceria is a promising catalyst for CO2 
activation and further reaction. 

IV. CONCLUSIONS 

We have performed analyses of structural geometries, energetics and electronic properties 
for the process of carbon dioxide adsorption to stoichiometric and reduced ceria surfaces, 
using the framework of density functional theory and employing the Hubbard U correction. 
It is found that oxygen vacancy plays an important role in the adsorption of CO2 on the 
ceria surface. While CO2 adsorption to stoichiometric Ce02 (HO) is relatively weak, with no 
significant changes to the geometry or electronic structure of the adsorbate, CO2 adsorption 
to reduced Ce02 (HO) is favored, particularly in the case of parallel CO2 adsorption atop 
the vacancy oxygen site neighboring two reduced Ce 3+ ; the energy of adsorption is -1.223 
eV. This chemisorption results in a bent CO2 molecule with an O-C-0 angle of 136.9° and 
elongated C-0 bonds compared to the gas-phase molecule. Bader charge analysis confirms 
that a net charge of —0.955 |e| is transferred from the reduced Ce02 (HO) surface to the 
adsorbed CO2 molecule in the most stable adsorption structure. Thus, the adsorbed CO2 
is activated to form the unidentate carbonate anion with the net negative charge localized 
on the oxygen atoms of the molecule. The excess spin density shows that the reduced 
ceria surface is partially re-oxidized. The re-oxidized Ce 4+ atom is closer in distance to 
the nO oxygen atom that moved towards the vacancy site upon surface reduction, while 
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the Ce 3+ atom is further from this nO oxygen atom and the adsorbed carbonate anion due 
to electrostatic repulsion. We suggest that reduced Ce0 2 (HO) is a promising catalyst for 
either thermal or electrochemical reduction of CO2 as the first step in the formation of fuels 
and chemicals for the closing of the carbon cycle. 
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